https://ogcruz.github.io/Bookdown/dados-de-area-iii.html
Exemplo com os dados de dengue em Dourados/MS
Nesta apresentação serão utilizados os dados que fizeram parte da dissertação de Isis Rodrigues Reitman intitulada “SAÚDE E AMBIENTE URBANO: A RELAÇÃO DE INCIDÊNCIA DE DENGUE E AS DISPARIDADES ESPACIAIS EM DOURADOS – MS”, apresentada no Programa de Pós-Graduação em Geografia da Universidade Federal da Grande Douradosos/MS, em abril de 2016.
Essa aplicação também se encontra no Curso de Estudos Ecológicos ministrado para o curso de Pós-Graduação em Epidemiologia em Saúde Pública em 2019, pelos pesquisadores Oswaldo Gonçalves Cruz (PROCC/FIOCRUZ) e Wagner Tassinari (DEMAT/ICE/UFRRJ)
Biliotecas do R que serão utilizadas
library(readr)
library(tidyverse)
library(sf)
library(maptools)
library(spatstat)
library(tmap)Lendo a tabela da população por setor censitário e os shapes files do contorno e por setor censitário de Dourados/MS
pop2010 <- read_csv('dados/pop2010.csv')
setor.sf <- read_sf('dados/Setor_UTM_SIRGAS.shp', crs = 31981)
contorno.sf <- read_sf('dados/contorno.shp', crs = 31981)Fazendo um join com as tabelas com os setores censitários + popoulação
setor.sf <- setor.sf %>% mutate (idsetor = as.numeric(CD_GEOCODI)) %>% inner_join(pop2010,by='idsetor')Lendo e plotando os casos de dengue georreferenciados em Dourados/MS
casos <- read_csv('dados/dengue_dourados.csv')
casos.sf <- st_as_sf(casos, coords = c("X", "Y"), crs = 31981)ggplot(setor.sf) +
geom_sf(fill = 'white', color='black') +
geom_sf(data=casos.sf, color='red',size=1) +
theme_void()Lendo e plotando os os pontos de coleta de lixo georreferenciados em Dourados/MS
lixo <- read_csv2('dados/lixo_dourados.csv')
lixo.sf <- st_as_sf(lixo, coords = c("X", "Y"), crs = 31981)ggplot(setor.sf) +
geom_sf(fill = 'white', color = 'black') +
geom_sf(data=lixo.sf,color='blue',size=1) +
theme_void()Como podemos observar, existem alguns pontos de coleta de lixo fora do contorno de Dourados/MS
Uma forma de ficarmos só com os pontos dentro do polígono é utilizando o comando st_intersection.
lixo2.sf <- st_intersection(contorno.sf, lixo.sf)ggplot(setor.sf) +
geom_sf(fill = 'white', color = 'black') +
geom_sf(data=lixo2.sf,color='blue',size=1) +
theme_void()Verificando as paletas de cores
library(RColorBrewer)
display.brewer.all()Utilizando as informações dos casos (pontos) + do lixo (ponto) + população de cada setor censitário (mapa temático)
ggplot(setor.sf) +
geom_sf(aes(fill=pop)) +
geom_sf(data=casos.sf,color='red',size=0.7, aes(colour = "Caso"),
show.legend = "point") +
geom_sf(data=lixo2.sf,color='salmon',size=1, aes(colour = "Lixo"),
show.legend = "point") +
scale_fill_distiller(palette ="PuBu", direction = 1) +
scale_colour_manual(values = c("Caso" = "red", "Lixo" = "slamon")) +
theme_minimal()Agora serão construidos buffers de 500m de distância ao redor de cada ponto de coleta de lixo. Isso é interessante para verificar se os casos de dengue ocorrem em um raio de até 500m de distância dos pontos de coleta de lixo. Ou seja, a pergunta é, será que a distância dos pontos de coeta de lixo influenciam a ocorrência do caso de dengue ?
Buffers: São polígonos que contornam um objeto a uma determinada distância. Sua principal função é materializar os conceitos de “perto” e “longe”.
lixo_buffer <- st_buffer(lixo2.sf, 500)
ggplot(setor.sf) +
geom_sf(aes(fill=pop)) +
geom_sf(data=lixo_buffer,color='gray', fill = "transparent", size=0.4) +
geom_sf(data=casos.sf, color='red',size=0.7) +
scale_fill_distiller(palette ="PuBu", direction = 1) +
scale_colour_manual(values = c("Caso" = "red", "Lixo" = "slamon")) +
theme_minimal()Represntando os casos e o lixo de forma interativa
tm_shape(setor.sf) + tm_borders("black") +
tm_shape(casos.sf) + tm_dots("red") +
tm_shape(lixo.sf) + tm_dots("green") +
tm_shape(lixo_buffer) + tm_borders("blue") +
tmap_mode("view")